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We formulate a model of the time-resolved Kerr rotation experiment on an ensemble of inde- 
pendent holes in a semiconductor nanostructure (e.g., confined in a quantum dot or trapped in a 
quantum well) in a tilted magnetic field. We use a generic Markovian description of the hole and 
trion dephasing and focus on the interpretation of the time-resolved signal in terms of the micro- 
scopic evolution of the spin polarization. We show that the signal in an off-plane field contains 
components that reveal both the spin relaxation rate and the spin coherence dephasing rate. We 
derive analytical formulas for the hole spin polarization, which may be used to extract the two 
relevant rates by fitting to the measurement data. 
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I. INTRODUCTION 

In recent years, considerable experimental progress 
has been achieved in the optical control and readout 
of spin states of electrons and holes in semiconductor 
nanostructures 1 - - — . Extended life times of spin states in 
quantum dots (QDs)2 and quantum wells (QWs) 1 ^ seem 
very promising for applications in spintronics, e.g. in 
the form of spin memor y 10 ' 11 , or in semiconductor spin- 
based quantum computin g 12 ' 13 . Particular expectations 
are related to hole spins, since the reduced hyperfine in- 
teraction in this case makes the hole spin decoherence 
even slowe r 14 ' 15 . 

In order to design feasible spin-based devices, one ob- 
viously has to understand the properties of spins in con- 
fined semiconductor systems. The most essential param- 
eters that must be known in order to predict the evolution 
of a spin in a real structure are the Lande g-tensor, which 
defines the unitary evolution of a spin in a magnetic field, 
and the relaxation (dephasing) times, accounting for the 
impact of the environment. As the Zeeman shifts can 
be rather small and, therefore, hard to resolve spectrally, 
many experiments rely on time-resolved methods 1 ^— . 
Among those, time-resolved Faraday 5 * 2 ^— or Kerr— ~— 
rotation measurements have proven to be very useful. In 
these experiments, one studies the rotation of the polar- 
ization plane of transmitted or reflected radiation (probe 
pulse) due to spin polarization excited by a circularly 
polarized pump pulse. In this way, the decay of spin 
polarization and the spin precession around the quanti- 
zation axis can be observed as a function of delay time 
between the pulses. 

The theoretical challenge in the description of confined 
hole spin properties is essentially three-fold. First, one 
has to model the nanostructure in order to find the effec- 
tive (/-factor for holes in a two- or zero-dimensional con- 
finement. This is usually done using k -p method a 18 ' 28 ' 29 
and the results are in reasonable agreement with mea- 
surements. Second, one needs to describe the spin de- 
phasing and relaxation processes. Phonon-mediated pro- 
cesses are often invoked for QD o 30 ' 31 , although other 



channels are also taken into account 3 ^. In QWs, phonons 
are expected to dominate for trapped holes, while scatter- 
ing due to compositional disorder is invoked for delocal- 
ized ones 3 - 3 . Another reason for dephasing may be system 
inhomogeneity, in particular g-factor fluctuations 34 . 

In the present paper, we deal with the third aspect of 
the problem, namely, the microscopic origin of the mea- 
sured signal and, more importantly, the relation between 
the magnetic (spin) orientation, which is supposed to be 
studied, and the optical field, which is experimentally 
accessible. We analyze also in what way the optical re- 
sponse (in particular, Kerr rotation) of the system de- 
pends on the parameters characterizing the system evo- 
lution. In the case of a time-resolved experiment with 
pulsed excitation, the physical interpretation of the de- 
tected optical response becomes non-trivial and consti- 
tutes a subject of study by itself 35,36 . This kind of the- 
oretical discussion has been presented for excitons in a 
QW— and, on a phenomenological level, for an n-doped 
QW system 2 ^. Very recently, a complete analysis of the 
Kerr and Faraday response for ensembles of n-doped QDs 
in an in-plane magnetic field was presented 3 - 6 -. Here, we 
focus on trapped hole states in QWs and on holes con- 
fined in QDs. We discuss the microscopic origin of the 
time-resolved Kerr rotation (TRKR) signal in a pump- 
probe experiment on a p-doped sample in tilted magnetic 
field. We perform a complete analysis in the density ma- 
trix formalism and describe hole spin dephasing on a gen- 
eral level, assuming only its Markovian character. Our 
description is applicable to various decoherence processes 
that have a well-defined Markov limit wchich is applica- 
ble under given conditions. Examples of such processes 
include phonon-assisted transitions or Coulomb scatter- 
ing. We discuss how the "longitudinal" and "transverse" 
dephasing rates (defined with respect to the tilted quan- 
tization axis) manifest themselves in the detected TRKR 
signal. 

As a result of our study, we point out that the optical 
signal follows the spin polarization in the limit of strong 
dephasing of optical coherences. The same holds true 
in the case of long-lived optical coherences if the phase 
relation between the pulses can be considered random 
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or, for an extended system, if the experimental geometry 
assures that the coherent part of the signal is emitted 
in a different direction. We show that the dephasing of 
the hole spin precession beats is governed by two rates 
involving a combination of the two relaxation rates. In a 
homogeneous system, this may allow one to extract both 
rates from a single experiment. In addition, we study 
the effect of inhomogeneity of ^-factors on the observed 
TRKR signal. 

The paper is organized as follows. In Sec. QJ we de- 
scribe the system and the experiment to be modelled. 
The next Sec. Mil contains the microscopic derivation of 
the optical TRKR signal. In Sec. IIV| we study the spin 
dynamics using a general model of Markovian decoher- 
ence, including also inhomogeneous dephasing. The rela- 
tion between the optical signal and the spin polarization 
is discussed in Sec. [V] In Sec. I VII we present the re- 
sults of our simulations and discuss the dependence of the 
TRKR signal on the off-plane tilt angle of the magnetic 
field and on the relative strength of different contribu- 
tions to dephasing (longitudinal, transverse, inhomoge- 
neous). These results are discussed in the next Sec. IVIII 
The appendix contains the derivation of the general Lind- 
blad equation describing Markovian decoherence in the 
hole-trion spin system. 



are the vectors of Pauli matrices corresponding to the 
hole and trion spin, respectively (the hole is treated as a 
pseudo-spin- 1/2 system). Here and in the following, we 
describe the system in a reference frame rotating with 
the zero- field hole-trion transition frequency uj. 

The spin states of each hole or trion can be described 
in terms of the "spin up" and "spin-down" states, that is, 
the basis states with definite projections on the growth 
axis (normal to the QW or to the plane of QDs), | f 
), 1 1) (for a hole) and |Tf>, \T I) (for a trion). For the 
magnetically isotropic trion, we define the two Zeeman 
eigenstates 



|T+) = cos-|Tt)+sin-m>, 
|T-) = -sin^|Tt)+cos^|T|). 

In the case of the hole, the quantization axis does not 
necessarily coincide with the field orientation. The hole 
spin eigenstates can be written as 



|+)=cos||t>+sin|||), 
-> = -sm^|t)+cos^>, 



II. THE SYSTEM 

We consider the optical response of a system composed 
of trapped (localized) holes which may be considered in- 
dependent (non-interacting). This situation may corre- 
spond to a p-doped quantum well in which holes are lo- 
calized by some kind of trapping potentials, e.g., inter- 
face fluctuations or nearby defects. The theory applies 
also to ensembles of quantum dots in a remotely p-doped 
structure. In thermal equilibrium, each trapping center 
is assumed to accommodate one hole. The area density 
of such trapped holes is v. Only heavy hole states are 
considered, since the heavy-light hole splitting is usually 
large in confined systems. The fundamental optical tran- 
sition at each trapping center consists in an excitation 
of an electron-hole pair which, together with the resident 
hole, forms a bound trion. It is assumed that the temper- 
ature is low and the driving pulses are spectrally narrow 
enough to restrict the description to the lowest hole and 
trion states. The system is placed in a magnetic field B 
oriented at an angle 9 with respect to the growth axis. 
The quantum well or layer of QDs is covered by a capping 
layer of thickness D. 

A single trapped hole-trion system is described by the 
Hamiltonian 

H = -^ B Bg h (T h - ^g t fi B B ■ er t , (1) 

where fi B is the Bohr magneton, g^ is the hole Lande 
tensor, g t is the Lande factor of the trion (i.e., essentially, 
of the electron), which we assume isotropic, and cr^cr^ 



where (j) is a certain angle depending on the structure of 
the hole Lande tensor. In our simulations, the latter will 
be assumed isotropic in the structure plane, so that 



at an 



<?h|| 



tan# 



(2) 



where g^± and <?h| are the in-plane and axial components 
of g. The Zeeman energy splitting for the hole is then 



huj^ = J g± sin 9 + gj| cos 2 0/ib-B, Note, however, that 

the actual structure of the Lande tensor enters the theory 
only via the the angle 4> and the Zeeman energy fkoh 
which can easily be found also for structures with a more 
complicated form of the Lande tensor—. 

A laser tuned to the trion line excites the system by 
two pulses (pump-probe configuration). The pulses prop- 
agate nearly perpendicular to the structure plane (a small 
deviation of the probe beam from the perpendicular axis 
is needed to separate the contributions to the third order 
response, as discussed in Sec.|V|. The first pulse arrives 
at t = and is circularly polarized (c+). The second, 
linearly (X) polarized one arrives at t = r. The ampli- 
tudes of the electric field in the two pulses (outside the 
semiconductor) are Ei = \Ei\e~ l ^ i , i = 1,2. The elec- 
tric field couples to the interband transitions via a dipole 
moment d. The pulse shape is described by a function 
7y(s), which is of the order of unity. The pulse length 
will be denoted by r p . The reflection amplitude at the 
semiconductor- vacuum interface is r = (1 — n)/(l + n), 
where n is the refractive index of the capping layer. The 
relevant Hamiltonian in the rotating wave approximation 
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is then 



+ ^=V (—) (ItXTtl + |I><T||) 



2v^ 
-H.c, 



(3) 



where Ai — d\Ei\{l + r) are the effective amplitudes of 
the two pulses inside the material. For the pulse shapes 
we will assume Gaussians, rj(s) = exp[— (l/2)s 2 ]. 

In addition to the evolution governed by the Hamil- 
tonian H — Hq + , the system undergoes dissipative 
dynamics due to the interaction with its environment. As 
long as this decoherence can be described in the Markov 
limit (which is reasonable in view of the relatively long 
time scales involved), its effect on the hole spin can be 
described by the universal Lindblad superoperator (see 
Appendix [A"]) . In this limit, the open system evolution is 
described by three dephasing rates: k± describe the lon- 
gitudinal decoherence, that is, spin relaxation between 
the Zeeman eigenstates for a given field orientation, while 
«o accounts for additional pure dephasing processes. The 
two spin-flip rates k± are related by 



k+ exp 



which guarantees the detailed balance condition at equi- 
librium. Here k + is the transition rate for a "down-flip" 
(from the upper to the lower Zeeman state), and k_ is 
the rate for an "up- flip". Note that, apart from this rela- 
tion, in a specific model of hole-reservoir interaction, the 
rates k± will depend on the Zeeman splitting and tem- 
perature via the spectral density of the relevant reservoir 
[see Eq. (|A6I) ]. However, in our general discussion, we 
treat one of them (or their sum) as a parameter of the 
model. 

The Lindblad dissipator describing the hole decoher- 
ence has the form 



Lh[p] = 



+ K- 



1 



(h) (h) i r (h) (h) \ 



. (4) 



where 



An analogous dissipator describes the spin dephasing of 
the trion. However, since the trion spin coherence time is 
much longer than its lifetime, the decay of spin coherence 
will be governed by the latter and the trion spin dephas- 
ing can be neglected. Note that the spin dephasing in 
the Markov limit is necessarily described in the eigen ba- 
sis |+), |— ) defined by the field orientation, which differs 



from the "spin-up" and "spin-down" basis defined by the 
structure symmetry and by the optical selection rules. 

The last part of the model is the radiative decay of the 
trion, which is accounted for by the Lindblad superoper- 
ator 



L Ta d[p] 



7i 



(t) (t) 
<7_ pa V 



(40 (4) 
-<r_ poT 



"2 70 



(5) 



where 71 is the radiative decay rate, 70 is the additional 
pure dephasing rate, and the transition operators are 



r (t) 



„(t) 



^0 



= ItXTtl 



Ai) 



r (4) 



i+xm 



ir4xm-itxti-ux+i 



Note that the distinction between the trion recombina- 
tion and pure dephasing is essential here not only because 
of the presence of various pure dephasing mechanisms in 
real systems^ - — but, much more importantly, because 
of the different effect these processes have on the spin- 
dependent optical response: both of them contribute to 
the decay of the optical polarization but pure dephas- 
ing, contrary to recombination, does not affect the trion 
spin occupations. In an ensemble of emitters, the de- 
phasing of trion coherences can be in fact dominated by 
inhomogeneous effects (distribution of the trion transi- 
tion frequencies). This would result in a different form 
of the coherence decay. However, from the point of view 
of the present study, this difference is of minor impor- 
tance and only the characteristic time of the coherence 
decay is essential. Therefore, we simplify the discussion 
by neglecting this kind of inhomogeneity and using only 
the pure dephasing rate 70 to characterize the optical 
dephasing. 



III. THE TRKR RESPONSE 

In this Section, we define the measured TRKR signal 
and clarify its relation to the microscopic variables (ele- 
ments of the density matrix) defining the state of the car- 
riers in a nanostructure at the moment when the probe 
pulse arrives. We show how the phenomenology of Kerr 
rotation emerges in the homodyne detection process from 
the interference of the macroscopic optical field reflected 
from the system surface with the radiation due to the 
interband optical polarization in the nanostructure. Fi- 
nally, we relate the latter to the spin polarization. 

The experimentally measured effect is a rotation of the 
polarization plane of the probe beam reflected from the 
sample. The total field is projected onto the two axes 
x, y, oriented at 45° with respect to the polarization of 
the probe beam. The rotation of the polarization axis 
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is given by the difference of intensity between the corre- 
sponding two components of the fiel d 24 ' 36 , 



only terms of the first order in the nanostructure response 
Eg± one finds the TRKR signal 



AI = — [<£*(*))-<£*(*))] = — Im (25+25-) > (6) 

where E + and 2?_ denote the (complex) amplitudes of 
the circularly right- and left- polarized components of 
the total field and (•) denotes time averaging over the 
period of the electromagnetic field. 

On the microscopic level, the observed reflected field is 
a sum of the beam reflected at the surface of the capping 
layer (this process will be treated on the usual, macro- 
scopic level) and the field emitted by the nanostructure. 
Thus, the two circular polarization components of the 
total field incident at the detector are 



E ± =E ] 



R± 



E 



S±j 



(7) 



where £r± is the field reflected from the surface of the 
capping layer and E$± is the field emitted by the carries 
trapped in the nanostructures. For pulsed excitation, 
slow evolution of the field amplitudes has to be taken 
into account. The field reflected at the surface simply 
follows the pulse envelope and the amplitudes of its cr+ 
and <t_ components at the sample surface are both equal 
to 



E K± {t) = E K {t) = ^j=rE 2V 



(8) 



The optical signal emitted from the structure origi- 
nates from the interband polarization. If p(t) denotes the 
density matrix representing the system state then each 
trapped hole-trion superposition contributes a cr+ com- 
ponent of the dipole moment 7r + = d(T^\p(t)\^)e~ iu)t + 
c.c. and a er_ component 7r_ = d(Tl \p(t)\ l)e~ lut +c.c. 
This results in the polarization current 

J+ = -ivud(T-\\p{t)\^)e-' luJt +c.c. 

(and analogous for J71 ) and the amplitudes of the corre- 
sponding two components of the radiation emitted from 
the structure (at the sample surface) are 



Es+(t) 
E S -(t) 



(9) 



where cp = 2Dujn/c is the phase shift (with respect to the 
field 2Sr) due to propagation through the capping layer. 

Note that if the two components have equal phases (as 
is indeed the case, see below) then the radiation emitted 
from the nanostructure is, in general, elliptical but its 
polarization axis is not rotated. Moreover, the intensity 
of this signal is weak. What one really measures in the 
homodyne detection scheme is, however, the signal co- 
herently superposed on the much stronger field reflected 
from the surface of the sample. Substituting Eqs. © 
and © into Eq. ([7]) and then into Eq. ^ and retaining 



AI(t) 



1 

p c 
1 



lm[E* R (t)E s+ (t) + E R (t)E* s _(t)] 



V ojdRe[E^(t)(T^\p(t)\^e-^ 
-E R (t)(Ti\p(t)\i)*e^] 



(10) 



In the above discussion, we have assumed that all the 
hole-trion systems evolve under the same conditions. The 
effect of inhomogeneity will be treated in Sec. IIVBI 
Eq. (|10l) describes the measured signal in terms of the 
quantum state of a nanostructure. This equation can 
be used to find the system response without any further 
simplifying assumptions based on a numerical simulation 
of the open system evolution. 

The next step is to derive the relation between the el- 
ements of the density matrix and the spin polarization 
before the arrival of the probe pulse. This relation can 
be expressed in an analytical form 36 under the assump- 
tion that the probe pulse is much shorter than any rele- 
vant time scale of the system dynamics (consistent with 
the idea that it is supposed to probe the instantaneous 
state of the system). One has to assume also that the 
dephasing times of interband coherences are longer than 
the pulse duration. 

In order to relate the Kerr response to the density ma- 
trix formalism we note that the system state p(i), which 
gives rise to the measured polarization, is prepared by the 
probe pulse from the state just before this pulse, p(r~), 
where t~ denotes the time instant just before the arrival 
of the probe pulse. Under conditions stated above, we 
can completely neglect the system evolution during the 
pulse. Then, the system density matrix is transformed 
according to p(t) = W (t) p(r~ ")W'(t) , with the unitary 
operator 

W {t) = cos 

-i sin ^ [(I t )(Tt | + 1 1) (Tl |) e*** + H.c] , 



where 



With this time evolution operator one finds for the 
interband matrix elements (for a =t, \) 



(Ta\p(t)\a) = 

„„„2 *2(t) 



(Ta\p(r-)\a) 



+ -sm^ 2 (t) [(Ta\p(T-)\Ta) (a\p(r-)W)] e"* 



. s ^^{a\p(r-)\Ta)e-^ 



(11) 
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Let us first assume that the delay time between the pump 
and the probe pulse is much longer than the interband 
dephasing time. (We will come back to the case when 
this is not fulfilled in Sec. [Vj) . In this case, the interband 
matrix elements at time t~ are negligible and we only 
have the contributions proportional to sin$2(i), i.e., the 
terms involving the occupation differences between trion 
and hole states. Substituting this into Eq. (|10[) and using 
the expression (j8]) for the reflected field we find for the 
TRKR signal 



AI 



1 



4V2 
x [St( 



rE 2 r] 



vuid sin ip sin $2 (t) 



Sh(r-)] 



(12) 



where 



Et(t) = (Tt\p(t)\Tt)-(Tl\p(t)\Tl), 

Sh(t) = (t|p(*)lt>-Ulp(*)U>. 



(13a) 
(13b) 



are trion and hole spin polarizations, respectively. 

We neglect here the delay between the field envelopes 
due to propagation through the capping layer, which is of 
the order of 1 fs, that is, much shorter than the picosec- 
ond pulse duration. The signal described by Eq. (fT2|) is 
proportional to the difference of hole and trion spin po- 
larizations just before the probe pulse. In this way, the 
TRKR measurement gives access to the evolution of the 
spin polarizations in the system. 

For a pulsed excitation, AI depends on time. We de- 
fine the time-integrated (TI) TRKR signal as 



A/ti = 



dtAI{t). 



(14) 



This quantity is a function of the time delay r between 
the pump and probe pulses. Since the homodyne re- 
sponse is proportional to the envelope of the probe pulse, 
the integration in the above equation is done over the du- 
ration of the probe pulse. 
We note that 



1 t — t \ \J~2h . 

dtrj [ — — — ] sin$2(i) = — ; — [1 — cos a.i\ 



where 



-4, 



-00 V T p 



V2h 



is the area of the probe pulse. Hence, the integrated 
detection signal is (see also Ref. |H) 



TI 



r sin tp 
4(1 + r) 



(l-cosc*2) [S t (T-) - E h (r")] 



In the weak pulse limit, the spin polarizations are pro- 
portional to the intensity of the pump pulse, hence the 
signal is also proportional to this intensity. Moreover, it 



follows directly from Eq. (|T5|) that in this limit the re- 
sponse is also proportional to the intensity of the probe 
pulse. 

The quantity vfko sets the natural energy scale for the 
emitted radiation and is equal to the energy the system 
would emit per unit area if each hole-trion system gener- 
ated one photon. Thus, AIth / (vhjj) , which is the quan- 
tity to be plotted based on the results of simulations in 
Sees. |V] and IVI1 corresponds to the average number of 
photons per one hole-trion emitter and one repetition of 
the experiment. 



IV. HOLE AND TRION SPIN DECOHERENCE 

In this section we present a detailed analysis of the spin 
dynamics based on an analytical solution of the equation 
of motion for the density matrix in an idealized situation 
of coherent optical driving and fast dephasing of optical 
coherences. The validity of these assumptions and the 
relation between the quantities calculated here and the 
actual signal are discussed in Sec. [V] We begin with 
a discussion of dephasing of a single system. Then we 
study the effect of inhomogeneity of g-factors across the 
ensemble. 



A. Homogeneous dephasing 

The set of dynamical variables describing the evolution 
of the system consists of the trion and hole populations 

N t (t) = (Tt|p(t)|Tt>+(T±|p(t)|T4,), N h (t) = 1-N t (t), 



the trion and hole spin polarizations defined in Eqs. (|13a[) 
and (|13bl) . as well as trion and hole spin coherences, 

X t (t) = (Tt\p(t)\Ti) + (Ti\p(t)\Tt), 
Y t (t) = i({Tt\p(t)\Tl)-{Tl\p(t)\Tt)), 
X h (t) = (t|p(*)U) + aiP(*)lt>, 

Y h ( t ) = i({t\p(t)\i) - a\p(t)\m 

Initially, all trion variables are zero. The spin of the 
trapped hole is in the thermal equilibrium state, which, 
in the basis of the hole spin eigenstates, is characterized 
by a spin polarization 



V = (+|p cq |+> - HPcqh) = tanh 



2k B T 



where p cq is the density matrix for the system state at 
equilibrium. This corresponds to the following initial val- 
ues for the dynamical variables of the holes: 

Nh\t<o = 1, Sh|t<o =pcos0, 
X h \ t<Q = psin<p, Yh|t<o = 0. 

In the present discussion, we assume that the pulse 
durations are much shorter than any characteristic time 



6 



scale of the spin dynamics and their action may be ap- 
proximately treated as instantaneous. This means that 
the excitation is coherent (the case of dephasing times 
comparable with pulse durations is discussed in Sec. fV|). 
We assume also that it is resonant (effects of detuning for 
the case of n-doped structures and exact Voigt geometry 
have been studied in Ref. [36h . Then, the effect of the 
pump pulse is to perform the rotation p — > V pV\ where 



V 



and 



11X11 



mxm+cos^dtxti + iTtxm) 



,Mn T (|tXTt|e* + |TtXt|e-*), 



Oil 
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FIG. 1: (Color online) The amplitudes of the contributions 
to the TRKR signal [Eqs. {17]) and {T8}] as a function of the 
tilt angle of the magnetic field. The amplitudes Ci and C2 
include both the hole contributions [Eq. (|18l) ] and the cor- 
responding trion contributions from Eq. (|17[) with the same 
time dependence. 



is the pulse area. The pulse generates the trion popula- 
tion and depletes the hole population accordingly, 



JVt(O) 



2 Oil pCOS( 



1 



2 2 

generates the hole and trion polarization 



N h (0) = 1 - JVt(O), 



Sh(0) 
St(0) 



(1 



cos 



2 ai 



) p COS <f> 



cos — sin 



2 oi 



oti pcos<p + 1 



and reduces the hole spin coherence which exists at ther- 
mal equilibrium in a tilted magnetic field 

Oil 

.Xh(O) = cos— p sin 0. 

The other dynamical variables remain zero. 

The subsequent dynamics of the system is generated 
by the Zeeman Hamiltonian Hq and by the dissipators 
L h and L rad , 



For the hole variables, the equation of motion ([16 
leads to the system of equations with iVh and E t actinj 
as source terms, 



1 + cos 



kq sin <f> + K\ - 



~z (kq - -y ) sin 2< t>Xh - Wh sin cf>Y h 
-k' cos (f>N h + 71 S t , 



X\i = — I k cos 



- — cos 2(f) j Xh + ujh cos 4>Yh 



11 
2 



Y"h — cos 4>Xh + ujh sin 0Eh , 



where k± = 



and 



P = -^[H Q ,p] + L h [p] + L T [p]. 



(16) 



For the occupation, this yields a single decay equation 
with the obvious solution 

N t =N t (0)e-^\ N h = l-N t . 

When the trion spin dephasing is neglected the trion vari- 
ables evolve according to a closed set of three equations, 

S t = -71 S t - Lo t sin6>Yt 
X t = -71 X t + uj t cos 0Y t 
Y t = -71 F t - cj t cos 6X t + uj t sin #S t , 

where u> t — gtp^B /h is the trion Larmor frequency. The 
solution for the trion spin polarization is easily found to 
be 



S t (i) = S t (0)e" 7lt [cos 2 6> + cosu;tism 2 0] 



(17) 



Eh = Eh — p cos cf>, Xh = Xh — p sin ( 



(the equilibrium values are subtracted). By the Laplace 
transform technique, one finds the solution for the hole 
spin polarization in the form 



(18) 



k=l 



The exponents can be found exactly in a simple form, 
while the expressions for the amplitudes Ck are rather 
lengthy, therefore we use the fact that the hole dephasing 
rates k±,o are much smaller than all the other rates and 
frequencies and give the formulas for Ck in the leading 
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Electron g factor 


9e = 


0.26 


Hole g factor 






- axial 


311 = 


0.6 


- in plane 


gx = 


= 0.04 


Trion recombination time 


1/71 


= 50 ps 


Refractive index 


n — 


3.44 


Pulse duration (pump & probe} 


r P = 


1 ps 


Pulse amplitude 






- pump 


d£i 


= 0.5 meV 


- probe 


dE 2 


= 0.1 meV 


Temperature 


T = 


1.6 K 


Magnetic field 


B = 


7 T 



TABLE I: System parameters which are fixed throughout the 
paper. The hole g factors are taken as for a 6 nm thick 
quantum wel l 18 ' 43 . Parameters correspond to an AlGaAs 
structure^. 



order, that is, neglecting terms 0(k±a). The result is 
Ai = -71, 
Ci 
A 2 



C 2 = 



2 w£ + 7i 
-iut - 71, 
1 71 wgcos 2 </>+ (71 +iu> t ) 2 



2 71 + iu t cj 2 + (71 + iu t ) 2 
A3 = 



sin 2 0E t (O); 



ft 



1 ui 



1 



■ cos 2 </>sin 2 0E t (O) + - sin 2^X^(0); 



2^-7? T 4 



«1 



A 4 = -iuj h - k - — , 

2 71 - iwh w t 2 + (71 - iw h ) 2 



G 



1 71 w t 2 cos 2 + (71 - iuifrf 



sin 2 </>E t (0) 



-sin 2 0E h (O)-~sin20X h (O); 



A 5 = 0, 

C5 = pcoscj). 

The values of the amplitudes Cj as a function of the ori- 
entation of the magnetic field for the parameters assumed 
in this paper (Tab. QJ are plotted in Fig. [TJ 

According to Eqs. (fl7|) and (|18p. there are three kinds 
of contributions to the total spin polarization Eh — E t . 
The constant one, (ftjAs), corresponds to the equilib- 
rium spin polarization. Exponentially decaying contri- 
butions, given by the first term in Eq. (|17|) and by the 
1st and 3rd term in Eq. (|18p. originate from the decay of 
the spin population with respect to the respective quan- 
tization axes. Since we assumed that the trion spin life- 
time is limited by the recombination time, the trion spin 
population decays with the recombination rate 71. The 
hole population decays with the spin relaxation rate K\. 
Since only spin polarization along the growth axis is rel- 
evant in the optical measurement, these contributions 
vanish when the spin quantization axis is perpendicu- 



lar to the structure symmetry axis, that is, cos0 = 
and cos = for the hole and trion contributions, re- 
spectively. It should be noted that due to the strong 
anisotropy of the hole g factor, the out-of-plane compo- 
nent of the hole spin is large already in a slightly tilted 
magnetic field. Therefore, the occupations of the Zee- 
man levels and their thermalization affect the optical re- 
sponse already in slightly tilted fields. The second term 
in Eq. ([17]) and the 2nd and 4th terms in Eq. JTH} re- 
flect the spin precession around the quantization axis. 
This precession affects the optically detected spin polar- 
ization only if the quantization axis is tilted with respect 
to the structure axis, that is, sin</> ^ and sin0 7^ 0. 



B. Inhomogeneous effects 

If the measured signal originates form an ensemble of 
emitters, it becomes dephased due to variation of sys- 
tem parameters across the ensemble. In our case, non- 
uniformity of g-factors makes the individual spins precess 
with various rates, which destroys the overall spin polar- 
ization. 

We assume that the number of systems in the ensemble 
is sufficient to describe the distribution of g-factors by 
a continuous distribution function. We neglect possible 
variation of the quantization axis. It is convenient to 
describe the inhomogeneity in terms of the trion and hole 
precession frequencies o>t and Wh, for which we assume 
Gaussian distributions 



1 



TZ? — 



t,h, 



(19) 



where uj t , uJh now become the central frequencies of the 
corresponding distributions. We assume that -C Wj, 
so that a variation of the amplitudes C\ in Eq. (fT8| can 
be neglected. Then, upon averaging with the distribu- 
tion functions (|19[) . the trion and hole spin polarizations 
[Eqs. dTTJ) and ([T8]> ] become 



St(*) 

Eh(*) 



E t (0)e 



-71* 



+ e 



5 

E 

fc=i 



fk(t) 



+ c.c. 



cos uj t t sin 2 0(2,0) 



(21) 



where the amplitudes Ck are the same as in Eq. (I18p . 
f k [t) = X k t for k = 1, 3, 5, f 2 (t) = -iutt - 7it - a 2 t 2 /2, 
and fi(t) = —iu>ht — («o + «i/2)i — <r 2 t 2 /2. As usual, 
the exponential decay of a single system is replaced by 
a Gaussian one if the dispersion of frequencies is larger 
than the homogeneous dephasing rates. 



V. SPIN POLARIZATION AND TRKR 
RESPONSE 

The discussion presented in the previous sections was 
based on some simplifying assumptions. On the one 
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hand, in our discussion of the TRKR response in Sec. IIIIl 
we concentrated on delay times longer than the interband 
dephasing time. Therefore we did not discuss the con- 
tributions to the interband polarization resulting from 
the interband coherences created by the pump pulse. In 
quantum wells, such interband coherences vanish very 
quickly but in self-assembled quantum dots their lifetime 
may be limited only by the recombination time, which is 
of the order of a nanosecon d 44 ' 45 . On the other hand, the 
derivation of the analytical formulas in Sec. HVl as well 
as of the relation between the spin polarization and the 
TRKR signal, is based on the assumption of coherent ex- 
citation. This, in turn, requires the coherence time to be 
long enough and the coherence assumption breaks if de- 
phasing of the optical coherences is fast. In this section, 
we will deal with these issues. 

In order to model the full optical response of the system 
we will numerically solve the evolution equation 

P = - T i H o + #ias , p] + L h [p] + L rad [p], (22) 
n 

calculate the optical signal according to Eq. (fTQ)) . and 
integrate the result according to Eq. (TT1)) . Some of the 
parameters will be kept constant for all the results pre- 
sented in this paper. The values of these fixed parameters 
(roughly correspopnding to an AlGaAs QW system sim- 
ilar to that studied in Ref. ;2j|) are collected in Tab. |U 
The trion Larmor frequency is Lu t — 0.16 ps -1 . The 
pulse amplitudes chosen here correspond to the pulse ar- 
eas a\ = 0.127T and ct2 = 0.016-7T for the pump and probe 
pulse, respectively. These values assure that the optical 
excitation is well in the linear regime, so that varying the 
pulse areas leads only to uniform rescaling of the signal 
intensity proportionally to the pulse intensities, that is, 
to a\ and a\. 

First, we will discuss the additional contributions in 
the case of delay times shorter or of the same order as 
the interband dephasing time. In this case, also the in- 
terband terms at time r _ in Eq. contribute to the 
total interband matrix elements at time t and thus to the 
emitted radiation. However, we notice that the phase of 
the second pulse ip2 enters differently in the three terms. 
Keeping in mind that the interband coherences created 
by the first pulse carry a phase factor e - *^ 1 (in the case 
of (Tt\p(r-)\t)) and e^ 1 (in the case of (t \p(T~)\Tt)), 
the total phases are ipi for the first term, -02 for the sec- 
ond term, and 2?/>2 — tpi for the third term. Hence, only 
the second term holds a fixed phase relation with the re- 
flected beam and can produce a non- vanishing homodyne 
signal if the relative phase between the pulses is random. 
Moreover, the two exciting laser pulses are usually ap- 
plied to the sample at slightly different directions fci and 
k>2- Then, for an extended system (ensemble of emitters), 
also the emitted radiation originating from the three con- 
tributions has different directions. The first one is emit- 
ted in the direction k\ of the reflected pump pulse, the 
second one in the direction of the reflected probe pulse &2 
and the third one in the background-free reflected four- 
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FIG. 2: (Color online) Comparison between the analytical 
result for the spin polarization (red solid lines) and the sim- 
ulated signal (dashed blue lines) for 9 — 86°, 1/ki = 100 ps, 
70 = (slow optical dephasing). In (a), the hypothetical de- 
tection signal corresponding to a fixed phase relation between 
the pump and probe pulses is shown; in (b) , the simulated sig- 
nal has been averaged over the phase of the probe pulse. 



wave mixing direction 2&2 — k%. Thus we find that the 
interband coherences resulting from the pump pulse ex- 
citation indeed give rise to emitted signals, but they do 
not contribute to the TRKR signal [Eq. ([TO])]. The only 
exception would be the case of temporally overlapping 
pump and probe pulses, where the actions of pump and 
probe pulses cannot be treated separately, and the case 
of perfectly aligned, phase-locked pulses where the inter- 
ference between the interband coherence from the first 
pulse and the second pulse gives rise to strong coherent 
control oscillations. Both cases, however, are not relevant 
for the purpose of this paper. 

The numerical solution for the system evolution ac- 
cording to Eq. (|22p does not involve any simplifying as- 
sumptions, apart from the resonant excitation condition. 
However, in this model, the geometrical relations between 
various directions in which the signals are emitted are not 
taken into account. Moreover, in a single simulation run, 
the relative phase between the pulses is fixed. Therefore, 
the calculated optical response contains both the TRKR 
signal and the coherent components. A comparison of 
the simulated signal to the spin polarization in this case, 
shown in Fig. [^a) , reveals that the former not only dif- 
fers by orders of magnitude but also is uncorrelated to 
the latter. This is the case even for delay times a few 
times longer than the trion relaxation time since the co- 
herent contribution belongs to a lower order of the opti- 
cal response and is many orders of magnitude stronger in 
the weak excitation limit. The coherent artefacts can be 
eliminated from the simulation result by simply averag- 
ing the results obtained with opposite signs of the probe 
amplitude. If this is done, the simulated signal agrees 
with the analytical formulas, as expected [see Fig. |U[b)]. 
This confirms that the approximations made in the ana- 
lytical solution do not noticeably affect the result. 

In the opposite case of strong interband dephasing, the 
analytical formulas are no longer valid. In Fig. 0(a) we 
compare the analytical result (red solid line) with the 
simulated signal for two values of the dephasing rate 70 
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FIG. 3: (Color online) (a) Comparison between the analytical 
result for the spin polarization (red solid lines) and the de- 
tection signal (dashed blue lines) for 9 — 86°, l/ni = 100 ps, 
70 7^ 0, as shown (fast optical dephasing case), (b) The sim- 
ulated signal has been rescaled up by a factor of 16 to show 
that its shape exactly follows the evolution predicted by the 
analytical formulas. 



which describes additional pure dephasing of the opti- 
cal coherence (beyond that associated with the radiative 
decay the rate of which, 71 is fixed throughout the pa- 
per). As the dephasing time becomes comparable with 
the pulse duration, the signal is quenched due to the 
reduced efficiency of optical pumping and probing. We 
note, however, that this quenching is uniform, that is, it 
does not modify the shape of the pulse. This is clear from 
Fig-EIb), where the simulated response for 70 = 5 ps -1 
has been multiplied by a factor of 16. Upon this rescaling, 
the simulated signal matches the analytically calculated 
one almost exactly. 

Thus, we have established the relation between the evo- 
lution of spin polarization in the system and the form of 
the TRKR response. It turns out that both the sim- 
ulated (or measured) signal and the analytical formula 
can yield consistent, correct information on the spin evo- 
lution. One has to eliminate the coherent polarization 
contributions from the calculated optical response in the 
slow dephasing case and the analytical formulas uni- 
formly overestimate the signal in the case of fast optical 
dephasing. 



VI. RESULTS 

In this Section, we discuss the evolution of the spin po- 
larization, based on the analytical solution to the equa- 
tions of motion derived in Sec. IIVI In all the simulations 
presented below, we set 70 = (hence, the term "pure 
dephasing" will always refer to the pure dephasing of spin 
states, described by the parameter kq). 

In Fig. S{a) , we show the evolution of the spin polariza- 
tion for a certain set of parameters. The signal appears 
to be dominated by two oscillating components. As dis- 
cussed in Sec. IIV A[ the short-period one corresponds to 
the trion precession with the frequency uj t . This contri- 
bution is damped with the rate 71 due to the finite trion 
life time. The other oscillating contribution originates 
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FIG. 4: (Color online) (a) The TRKR signal for 9 = 87°, 
l//ti = 100 ps, kq = Uh = 0t = 0. (b) The real part of the 
Fourier transform of the signal (red solid line) and the contri- 
butions from the trion precession (green short-dashed line), 
hole precession (grey dotted line) and hole spin relaxation 
(blue long-dashed line). The contribution from the exciton 
spin relaxation is invisible on this scale. 



from the hole precession and is damped with the total 
hole spin dephasing rate Ki/2 + kq, reflecting the decay 
of the hole spin coherence (transverse dephasing). Two 
other contributions, which are less evident in the plot, 
have a non-oscillating character and reflect the spin re- 
laxation leading to thermalization between the Zeeman 
eigenstates (longitudinal decoherence) . The presence of 
these parts of the signal becomes easily visible in the form 
of a central line in the Fourier transform of the TRKR 
response, shown in Fig. 0Jb). We plot here also the indi- 
vidual contributions following from Eqs. (fTTl) and (IT8)) . 

The two most interesting aspects of the TRKR re- 
sponse are the dependence of the signal on the tilt angle 
between the magnetic field and structure plane and the 
effect of various dephasing types (spin relaxation, pure 
dephasing, inhomogeneous dephasing). In the following 
subsections, we start our analysis with the angle depen- 
dence and later proceed to the role of various dephasing 
contributions. 



A. Tilt angle dependence 

Due to the strong anisotropy of the hole g factor, the 
TRKR signal shows a very strong dependence on the an- 
gle at which the magnetic field is tilted off the system 
plane. The quantization axis of the hole spin is far from 
the plane even for small tilting angles. Therefore, the 
form of the signal changes strongly when 9 is varied in 
the range of a few degrees from the in-plane orientation. 

In Fig. [5ja) we show the TI-TRKR response for three 
different tilt angles 90°— 8 between the magnetic field and 
the structure plane. In these calculations, we keep the 
same value of K\ for all angles even though the Zeeman 
splitting changes, which is to some extent artificial. A 
correct dependence would follow from the detailed mod- 
eling of a specific decoherence channel, which is beyond 
the scope of the present general description. In all the 
cases shown in Fig.[5ja), one can clearly see a similar con- 
tribution from the trion Larmor precession. However, the 
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FIG. 5: (Color online) (a) The TRKR signal for the three tilt 
angles 90° — 9 as shown, 1 /ki = 100 ps and kq = Oh = ft = 0. 
(b) The real part of the Fourier transform of the three signals 
(line coding as in (a)). Only the positive frequency part of 
the symmetric spectrum is shown. 



FIG. 6: (Color online) Comparison of the TRKR response (a) 
and its Fourier transform (b) for different contributions from 
the pure dephasing as shown (the line coding is the same in 
both panels). The value of Ko is adjusted so that fti /2 + kq = 
1/100 ps -1 in all three cases. Here 90° — 6 = 3°, a t = <Jh = 0. 



hole contributions are different. The amplitude of the os- 
cillations decreases as the field is tilted more off the plane. 
At 90° — = 10°, the hole contribution is dominated by 
a monotonous decay superposed on the trion oscillations. 
The reason for this is clear: For the 1° tilt, the hole pro- 
cesses around an almost in-plane quantization axis (ori- 
ented at about 15° off-plane). Such a precession leads to 
a strong variation of the perpendicular component of the 
spin, while the thermalization of the spin occupations is 
associated mostly with the optically irrelevant decay of 
the in-plane component. On the contrary, according to 
Eq. at 90° — 9 = 10° the hole spin quantization axis 
is close to perpendicular (90° — 4> — 70°). The precession 
then takes place mostly in the plane, while the spin pop- 
ulation decay affects the perpendicular spin polarization 
and is visible in the experiment. 

This qualitative difference in the system evolution is 
visible even more clearly in Fig. [S^b), where we plot the 
real part of the Fourier transform of the TI-TRKR sig- 
nals shown in Fig. 03a). Three characteristic features 
are visible in this spectrum. Starting from the right, the 
broad one at co = 0.16 ps -1 corresponds to the trion 
precession. The orientation of the magnetic field does 
not affect the position of this feature because the trion 
(electron) g factor is isotropic. Moreover, for the narrow 
range of tilt angles considered here, the effect on the am- 
plitude of the trion oscillations is very small. Therefore, 
this feature is almost insensitive to the orientation of the 
field in the considered range. The second feature moves 
from lo = 0.025 ps" 1 at 90° - 9 = 1° to lo = 0.07 ps" 1 at 
90° — = 10° and looses its amplitude. It corresponds to 
the hole precession. The frequency shift is obviously due 
to the growing contribution of the large axial component 
of the hole g factor as the field is tilted off the plane. The 
decrease in amplitude corresponds to the fast reorienta- 
tion of the hole spin quantization axis, which leads to 
reduced contribution of the hole precession to the optical 
signal. The third feature is the central line, correspond- 
ing to the exponential decay of the hole spin population. 
As the magnetic field is oriented more off-plane, the con- 
tribution of this process to the spin polarization grows 



and this feature becomes stronger. 

B. Dephasing 

Another interesting feature observed in the simulation 
is a different dependence of the decay time of two hole- 
related components on the two hole decohcrence rates kq 
and K\. This is visible in Fig. [51 where we fix the pre- 
cession damping rate ki/2 + kq and change the relative 
contributions from the spin relaxation (ki) and the addi- 
tional pure dephasing (kq). In the time-resolved picture 
[Fig. [Bfa)], the differences are not particularly charac- 
teristic, except for the long exponential tail which de- 
velops as the spin relaxation becomes very slow. Much 
more pronounced differences can be noticed in the Fourier 
transform [Fig. [HJb)] . As the parameter modification af- 
fects only the hole dynamics, the trion feature remains 
unchanged. Moreover, since we fixed the total dephasing 
time of the hole precession, the feature at the hole Lar- 
mor frequency, uj — Wh, changes very little (only due to 
the change in the tails of the neighboring zero-frequency 
feature). On the contrary, the central line changes very 
strongly. As the lifetime of the spin population becomes 
longer, this line gets narrower, with the line area remain- 
ing constant. It seems, therefore, that the spectral com- 
ponents of the TRKR signal in a tilted magnetic field 
carry useful information on the relative strength of dif- 
ferent contributions to spin dephasing. 

Another source of damping of the observed spin preces- 
sion oscillations in the case of an ensemble measurement 
is the inhomogeneity of Larmor frequencies due to a vari- 
ation of g factors of the individual emitters in the ensem- 
ble. Obviously, only the precession-related contributions 
are sensitive to the inhomogeneity effects. Indeed, as fol- 
lows from Eqs. (f2"0"|) and (f2"Tj) inhomogeneity affects the 
oscillating contributions, while the exponentially decay- 
ing ones remain unaffected. Since the main concern in the 
present work is the hole spin evolution, we will restrict 
the discussion to the case of a t = 0. In Fig. we show 
the evolution of the TRKR signal and the correspond- 
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FIG. 7: (Color online) Comparison of the TRKR response 
(a) and its Fourier transform (b) for homogeneous and inho- 
mogeneous dephasing. Only the frequency range relevant to 
hole dynamics is plotted in (b). Red solid lines: homogeneous 
pure dephasing, Ko = 3/400 ps _1 , <7h = 0. Blue dashed lines: 
inhomogeneous dephasing, k,q = 0, Oh = 0.0147 ps" 1 . In both 
cases, 90° - 9 = 3°, Ki = 1/100 ps" 1 . 

ing spectrum for a fixed value of the hole spin relaxation 
rate Ki with the dephasing contribution dominated by 
the homogeneous pure dephasing (red solid lines) and by 
g factor inhomogeneity (blue dashed lines). The param- 
eters Ko and <7h are chosen such that the full width at 
half maximum of the damping envelope (in the time do- 
main) is the same in the two cases. Again, there is only 
a minor change in the time-resolved signal [Fig. [7(a)]. 
However, the shape of the spectral feature correspond- 
ing to the hole spin precession changes from Lorentzian 
to Gaussian. This change may be characteristic enough 
to discriminate the homogeneous vs. inhomogeneous de- 
phasing in a real measurement data. 



VII. DISCUSSION AND CONCLUSION 

We have developed a complete theory of the time- 
resolved Kerr rotation experiment for a system of trapped 
holes in tilted magnetic fields. The theory is applicable to 
quantum dots or weak trapping centers in quantum wells. 
In our approach, we adopted a general description of hole 
spin relaxation and dephasing in the Markov limit, based 
on the Lindblad equation for the open system dynam- 
ics. Spin dephasing is a rather slow process so that the 
Markov approximation should work well for this problem 
and our approach can be expected to cover a wide range 
of physical effects in a way which is independent of the ex- 
act microscopic mechanisms. One should note, however, 
that there are dephasing mechanisms that do not admit a 
Markov approximation of this kind. The most important 
example of this class is a spin- environment coupling via 
Heisenberg-like (spin-spin) interaction Hamiltonian. 

Our analysis shows that the hole spin dephasing con- 
sists actually of two processes the relative contribution 
of which depends on the tilt angle between the mag- 
netic field and the structure plane, with an important 
role played by the strong anisotropy of the hole g fac- 
tor. These two processes are relaxation between the 



Zeeman states (occupation thermalization) , which dom- 
inates the optical response when the quantization axis 
is close to perpendicular to the plane (aligned with the 
structure axis), and dephasing of coherences between the 
spin states, which contributes mostly when the quantiza- 
tion axis is close to the plane. It should be kept in mind 
that the hole spin quantization axis is always much closer 
to perpendicular than the magnetic field orientation. 

As both these dephasing contributions are marked in 
the optical signal for a slightly tilted field (a few degree) 
a single set of experimental data conveys, in principle, 
the full information on the spin dephasing. Extracting 
this information is not straightforward at least for three 
reasons. First, the Larmor frequencies are not much 
higher than dephasing rates and the spectral features re- 
lated to these two dynamical contributions overlap rather 
strongly. Second, the coefficients of Eq. (fT8)l are complex 
and the features are not purely Lorentzian. Third, in 
the case of an ensemble experiment, inhomogeneous de- 
phasing can dominate the intrinsic one. On the other 
hand, Eqs. (fTT)) and (fT5|) provide analytical formulas for 
the spin polarization. As we have shown in Sec. Mil this 
spin polarization is identical with the measured signal 
(up to uniform rescaling). Then, the formulas provided 
by our theory can be used to fit the experimental data 
with just a few parameters, which might allow one to ex- 
tract all the relevant decoherence rates. Moreover, in the 
present paper, we have used parameter values which cor- 
respond to a quantum well system, where the dephasing 
is rather strong. In quantum dots, where spin coherence 
times are much longer, the signal should show much more 
pronounced and separated features, which can make the 
analysis much easier. Finally, as shown in Fig. [I] the 
imaginary parts of the amplitudes are relatively small, 
so even a rough line width estimate based on the Fourier 
spectrum of the time-resolved signal could yield reason- 
able information on the decoherence rates. 
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Appendix A: The Lindblad equation for the spin 
dephasing 

In this Appendix, we derive the general Lindblad equa- 
tion which governs the dissipative evolution of the den- 
sity matrix of a trapped hole in the Markov limit (an 
analogous equation can be written for the trion spin). 

Any observable related to the two-level spin system 
can be written as a combination of Pauli matrices a± t o 
acting on the Hilbert space of hole spin states and written 
in the basis of Zeeman eigenstates for a given orientation 
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of the magnetic field. Therefore, we can write the general 
Hamiltonian for the system-reservoir interaction in the 
form 



flint = 



^2 <nRi, 



(Al) 



i=±,o 



where Ri are certain operators on the Hilbert space of 
the reservoir and i?+ = R_ . 

One starts with the exact equation of motion for the 
reduced density matrix p of the hole spin in the interac- 
tion picture 



dp{t) 

dt 



1 



dr Tr R [H int (t), [H int (r), g(t)}} , (A2) 



where g(t) is the density matrix of the total system, Trjj, 
denotes partial trace with respect to the reservoir degrees 
of freedom, and to is the initial time of the evolution. 

Let us denote the reservoir memory time by r mcm . The 
Markov approximation is based on three assumptions 46 : 
(1) the time t of interest is much longer than r mcm ; (2) 
the change of the system state (in the interaction pic- 
ture) is small over the time r mcm ; (3) the relaxation of 
the reservoir to its thermal equilibrium is fast compared 
to the rate with which it is excited by the system evo- 
lution, so that the total density matrix of the system 
can be written in a product form, with the reservoir at 
equilibrium. Eq. (IA2[) can be then approximated as 

dp(t) = 
dt 

-i J ds Ttr [flint (t), [flint (t ~ S), p(t) ® p R |A3) 

where pr is the thermal equilibrium density matrix of 
the reservoir. 

In the interaction picture, we denote the reservoir op- 
erators Ri by Ri (t) and write the hole spin Pauli matrices 
as 07(f) = aie~ l0Jl , where w_ = — lu+ = uj a and ujq = 0. 
We define the reservoir spectral densities 



Rij{uj) 



2irh 2 



dte iujt (Ri{t)R 



31 > 



(A4) 



where (Ri(t)Rj) = Ttr pRRi(t)Rj. With this definitions, 
transforming Eq. (|Alj) to the interaction picture and sub- 
stituting into Eq. (|A3|I we get 



POO 

+ {paia 3 - afpai) \ dse^'+^s 
Jo 



In the next step, we use the identity 



1 



dse ±lUs =TrS(Q)±iV- 

where V denotes the principal value. Moreover, we note 
that the terms with w; + luj oscillate quickly in time 
and do not contribute considerably to the evolution of 
the density matrix. We can thus write 



dp{t) 
dt 



= 2-K^~5l] Rlj (Ujj ) ( <7j poi -^{<Ji<J j7 p} 



:[h,p], 



(A5) 



where 



h = hyi/P I du^i^-atc 



0, {A,B} + = AB + BA 



5ij = 1 if and only if uii + ujj 
and [, 1 denotes the commutator. 



The second part of the right-hand side of Eq. (|A5I) . 
containing the commutator, is a correction to the uni- 
tary evolution due to environment-induced level shifts. 
These effects are very weak and amount only to a small 
renormalization of the (/-factor. We will, therefore, dis- 
regard this term. Of interest to us is the first term, de- 
scribing the dissipative impact of the environment. It is 
clear that, irrespective of the nature of the reservoir, the 
dephasing in the Markov limit, in a given experimental 
situation, is completely described by three rates, 



k± = 2nR ±T (uj h ), k = 47i\R O o(0). 



(A6) 



However, using Eq. (|A4[) . it can be shown that 
R-+(-uj) = e~^/( feBT )i? + _(u;), where fc B is the Boltz- 
mann constant and T is the temperature. Hence, the 
number of dephasing parameters reduces to two. These 
two dephasing rates are related to the longitudinal and 
transverse dephasing times T\ and Ti (with respect to 
the quantization axis) by the usual formulas 



Ti = 



T 2 = 



1 
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